### Do all regressions:

####### Cleaning working space
rm(list = ls(all.names = TRUE)) 
gc()

## Read data:
data <- read.csv("intermediary_outputs/base_with_connection_data.csv")

# Create auxiliary functions:
source("09a_supp_function_regressions_en.R", encoding = "UTF-8")

# Regressions 1 - Table 2
source("09b_regs_1_en.R", encoding = "UTF-8")

# Heterogeneity by School Level, Gender, Class composition, Grades (S6)
source("09c_sep_het_linear_en.R", encoding = "UTF-8")

